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We identify a non-Hermitian chiral random matrix theory that corresponds to two-color QCD 
at high density. We show that the partition function of the random matrix theory coincides with 
the partition function of the finite- volume effective theory at high density, and that the Leutwyler- 
Smilga-type spectral sum rules of the random matrix theory are identical to those derived from the 
effective theory. The microscopic Dirac spectrum of the theory is governed by the BCS gap, rather 
than the conventional chiral condensate. We also show that with a different choice of a parameter 
the random matrix theory yields the effective partition function at low density. 
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Random matrix theory (RMT) has found numerous 
applications within and outside of physics, e.g., in nu- 
clear physics [1], mesoscopic condensed matter physics 
[2], quantum chromodynamics (QCD) and QCD-like the- 
ories [3], quantum gravity [4], number theory [5], econo- 
physics [6], wireless communications [7], and many more. 
Most of these applications are described by Hermitian 
random matrix models. While the first paper on non- 
Hermitian RMT by Ginibre appeared already in 1965 
[ ], there was essentially no activity in this area for a 
long time. However, in the past 10 years or so numer- 
ous new applications have been found in many fields of 
science that can be described by non-Hermitian random 
matrix models [9]. Prominent examples are QCD and 
QCD-like theories at nonzero density which, depending 
on their anti-unitary symmetries, are described by chi- 
ral versions of the Ginibre ensembles [10]. So far, the 
applications of non-Hermitian RMT in QCD have been 
restricted to moderate densities, roughly up to the crit- 
ical density for the chiral phase transition. However, in 
recent years there was quite some activity in the study 
of QCD and QCD-like theories at very high densities, 
by means of effective theories and renormalization group 
methods, and a number of interesting new phenomena 
such as color superconductivity [11] and color-flavor lock- 
ing [ ] have been elucidated. The present paper is the 
first application of non-Hermitian RMT at high density 
and is expected to open up many more applications. We 
shall see that the situation is quite different compared to 
the one at low density. We will focus on the case of QCD 
with two colors which can in principle be tested in lattice 
simulations because the fermion sign problem is absent 
(for pairwise degenerate quark masses). 

Recently we have constructed the low-energy effective 
theory for high-density QCD with two colors and an even 
number Nf of fermions in the fundamental representa- 
tion [13]. At sufficiently large quark chemical potential // 
(fi 3> Agu(2)), chiral symmetry is broken spontaneously 
by a diquark condensate as 

SU(JV» L x SU(N f ) R x U(1) B x \J(l) A 
^Sp(N f ) L xSp(N f ) R , (1) 



and the degrees of freedom of the effective theory are the 
Nambu-Goldstone (NG) bosons associated with (1). 

Considering the theory in a finite volume, we have 
identified the corresponding microscopic regime (defined 
by the requirement that the partition function is domi- 
nated by the zero-momentum modes of the NG bosons) 
and derived a set of Leutwyler-Smilga-type spectral sum 
rules [ 1 4] obeyed by the complex eigenvalues of the Dirac 
operator. We discovered that at high density the dimen- 
sionful quantity governing the Dirac spectrum is not the 
usual chiral condensate but the BCS gap A associated 
with the condensation of quark-quark pairs. This is also 
true for three-color QCD [ ]. 

The sum rules we obtained earlier plus a universal- 
ity argument hinted at the existence of a corresponding 
RMT. Here, we propose that the non-Hermitian chiral 
RMT describing two-color QCD at large /j, for an even 
number of fundamental fermions is given by 

Z({m/}) = J dAdB e -Ntr(AA^+BB T ) 

where A and B are real N x N matrices, the integration 
measure is the flat Cartesian measure, the rrif are dimen- 
sionless quantities corresponding to the quark masses, 
and Nf is assumed to be even. At this stage physical 
scales have not been introduced yet. (The matching of 
random-matrix quantities to physical quantities is given 
in Eq. (9) below. There is no need to introduce a Gaus- 
sian width parameter for A and B, as this parameter can 
be absorbed in the m/.) Note that A and B are square 
matrices. In principle, the model (2) can be extended 
to rectangular matrices A and B [16], see Eq. (4) be- 
low, giving rise to topological zero modes. However, at 
high density the topological susceptibility is strongly sup- 
pressed [17] so that topological zero modes are irrelevant 
in the physical situation we are studying. 

We stress that our model (2) differs from earlier RMT 
approaches to QCD at fi ^ [18-20]. In these ap- 
proaches, the corresponding low-energy effective theory 



2 



is formulated in terms of NG bosons parametrizing the 
coset space SU(2iV/)/Sp(2iV/) of the symmetry breaking 
pattern by a chiral condensate at fi = [21], i.e., /x is con- 
sidered as a small perturbation. On the other hand, the 
low-energy effective theory corresponding to (2) is for- 
mulated in terms of NG bosons parametrizing the coset 
space SU{N f ) x SXJ(N f )/[Sp{N f ) x Sp(N f )] (see [13] and 
Eq. (8) below), i.e., it realizes the pattern of spontaneous 
chiral symmetry breaking induced by BCS-type diquark 
condensation at high density, with a vanishing chiral con- 
densate. To see this explicitly, we rewrite (2) in the chiral 
limit as 

Z({0})= JdAdet N f/^_° AT e"^ T 

dct N >/ 2 (_° T fj e - Nt * BBT . (3) 

Each of the two factors on the RHS of this equation 
corresponds to a chiral orthogonal ensemble at /i = 0, 
for which the symmetry breaking pattern is known to 
be SV(N f ) Sp(Nj) [22]. Thus, we can immediately 
conclude that Eq. (3) exhibits the symmetry breaking 
pattern SU(iV» x SU(JV/) -> Sp(N f ) x Sp(iV/), which 
agrees with Eq. (1). There are also NG bosons associated 
with the breaking of U(1)b (baryon) and U(1)a (axial). 
The former decouples completely and does not affect the 
quark mass dependence of the partition function [13]. 
The latter decouples in the chiral limit in which (3) was 
considered. At nonzero quark masses, see (2), we have a 
nontrivial U(1)a integral, see Eq. (8) below. 

The model (2) is not new. It corresponds to class 
2P in Magnea's classification of non-Hermitian ensem- 
bles, see Table 2 of [ ■]. What is new is the realization 
that Eq. (2) describes dense two-color QCD. Below we 
show explicitly that the effective theory resulting from 
Eq. (2) at large N exactly reproduces the finite-volume 
partition function of dense two-color QCD in the micro- 
scopic limit, obtained in [13] based on symmetry prin- 
ciples and weak-coupling calculations at large fi. This 
strongly suggests that Eq. (2) is the correct RMT for 
dense two-color QCD, paving the way towards a quan- 
titative understanding of the Dirac spectrum at large /i. 
A full proof of the equivalence, which would also estab- 
lish the equality of all spectral correlation functions, re- 
quires the study of the partially quenched theory, see [24] , 
which we leave to future work. However, we will also see 
that the Leutwyler-Smilga-type sum rules derived from 
Eq. (2) coincide with those obtained in [ L ] . These sum 
rules are moments of the spectral correlation functions 
and thus provide a nontrivial piece of evidence for the 
equivalence. 

Before going into details, let us list some basic prop- 

/ A 

erties of the model (2). First, the matrix D = I ^ _ 

is neither Hermitian nor anti-Hermitian. This implies 
that the eigenvalues of V are distributed over the com- 
plex plane. Second, T> preserves chiral symmetry (i.e., 
it anticommutes with 75). Third, it follows from chiral 



symmetry and from the fact that A and B are real that 
the eigenvalues of T> appear either in pairs (A, —A) with 
A € K U iR or in quartets (A, -A, A* , -A*) with A e C 
[ ] . All of the above features are in common with dense 
two-color QCD [13]. 

Let us also elucidate the relation between our ensem- 
ble (2) and related ensembles in the literature. The first 
two-matrix model for fj, 7^ 0, in which the matrix multi- 
plying fi is not a unit matrix but an independent random 
matrix, was introduced by Osborn [ ] . His original anal- 
ysis was restricted to Dyson index j3 = 2, corresponding 
to three-color QCD. Subsequently, extensions to = 4 
[ ] and j3 = 1 [16] were constructed. Two-color QCD 
corresponds to /3 = 1, and the two-matrix model, which 
was studied mathematically in [28] for Nf = 0, reads 



Z v (p>, {r 



dCdD e -™*(CC T +DD* 



Nf 



Y[ det 



/=i 



m f 



c 



fiD 1 



[ID 

TO/1 



(4) 



Here, C and D are N x (N + v) real matrices and /t 
is a parameter corresponding to the physical chemical 
potential fi. For v — and fi — 1, (4) reduces to (2). 

We shall now try to derive the effective theory corre- 
sponding to (4) at large N for < jl < 1, following a 
similar analysis at fi = [22]. The mass terms in (4) 
are generalized to a generic mass term Pr + mP^ 
as in the QCD Lagrangian, where dim(m) = Nf and 
Pr/l — 5(1 ±75)- We introduce Grassmann variables to 
write the determinants as exponentials, integrate out the 
Gaussian matrix elements, and then introduce auxiliary 
matrices for a Hubbard-Stratonovich transformation. Af- 
ter integrating out the Grassmann variables we obtain 

Z v (fi, m) = J dKdLdP e -SNtr(K^ + LL^ + 2PP^ 



X Pf 



x Pf J 



,v ( v/i + A 1 ^ -VT^Wp t - |i 

y/1 - Pl 2 P + §mt v 71 + f VK 

N+v ( ji + Wtf -^y^Wp*-\ 

v^A^pt + i m yTTA2£ 



m 



(5) 



where K, L and P are Nf x Nf complex-valued matrices, 
with K and L antisymmetric. Pf denotes the Pfaffian of 
the matrix. So far no approximation has been made. 

The so-called strong non-Hermiticity limit is defined by 
N ^> 1 with jl fixed. We now focus on the case relevant 
for us, i.e., /I = 1 and v = 0. In this case P drops out 
and we find 



Z (l,m) = JdKdL e -^tr(KK^ + LLi) 



X Pf 



mt K ) \m L 



, (6) 



where the normalization has been changed slightly. 

Below we assume that Nf is even. Since K and L are 
antisymmetric they can be brought to the standard form 
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K = UAU T with U G U(iV / )/[Sp(2)] A '// 2 and A a real 
antisymmetric matrix with Afc^+i = — Afe+i^ > and all 
other matrix elements zero (likewise for L). For ||m|| <C 1 
the integration over A can be estimated by a saddle-point 
approximation at m = 0, while the integration over U is 
soft and has to be done exactly at m ^ 0. Elementary 
calculation shows that the saddle point is given by A = 
1 N 



I/y/2 with I = 



1 



Since U enters Zq only in the 



form UIU T , the integration manifold can be enlarged to 
\J(N f ). Thus 

Z„(l,m)~ J dUdV Pf*( ^ V mf 1{//[/T 



m 



X Pf 



m' 



J dUdV exp[- 2NRctr(mUIU T m T V*IV 1 )] 



U(JV» 



(7) 



where we expanded the exponent to the first nontrivial 
order in m. Redefining the variables as U = U e l ^ e+ ^ 
and V = VS 6 -'^ with U, V G SU(iV» and < <p < n, 
and setting A = e 2lip , we obtain 



Z (l,m)~ j dA j dUdV 

U(l) SU(W f ) 

x exp [ - 2A^Retr(^ 2 mi7/{7 T m T V r */y t )] , (8) 

which exactly coincides with the finite-volume partition 
function at large fi [13, Eq. (4.2), Eq. (4.24)] if we identify 



N m 



4vr 2 



(9) 



where M is the dimensionful physical mass matrix and V4 
is the space-time volume. Therefore, our model (2) de- 
pends on the chemical potential fi only implicitly through 
A. (Note that the gap A depends on [i via the relation 
A ~ [leT 1 ! 9 [ )] and that A is related to the diquark 
condensate via the relation (qq) ~ /J 2 A/g [13, 30].) The 
absence of terms linear in m in the exponent of (8) is 
a consequence of the Z(2)j, x Z(2)# symmetry (flipping 
left-handed and right-handed quarks independently, i.e., 
q L -> ±q L and q R -t ±q R ) of the diquark pairing (qLqh) 
and (qRqR) , which distinguishes the model (2) from con- 
ventional chiral RMTs. The physical limits in which the 
model (2) describes two-color QCD are given by the in- 
equalities defining the microscopic domain (or lowest or- 
der of the e-regime) at high density [13], 



1 T 1 

T « i « 

A m n,n,v' 



(10) 



where L is the linear extent of the Euclidean box and 
m Yl jj , is the mass of the NG bosons at high density. 



In addition to the mass term, we could also add a di- 
quark source term cx jqq to the model (2). It can be 
shown that the partition function at N 1 then con- 
tains a term linear in j, resulting in a nonvanishing di- 
quark condensate in this model. The addition of the 
diquark source term makes the model somewhat more 
complicated and will be addressed in future work [31]. 

Since in the large-iV limit the partition function (2) 
has the same mass dependence as the partition function 
(8) of the effective theory, it follows immediately from the 
calculations performed in [13] that the Leutwyler-Smilga- 
type sum rules for inverse Dirac eigenvalues of the model 
(2) agree with those obtained in [13], as stated earlier. 

For Nf = 2 we can also calculate the partition function 
at finite N. In [10] the correlation function of two charac- 
teristic polynomials was studied in the quenched version 
of the model (4), which corresponds exactly to Z v for 
Nf = 2 in (4). The result for fi = 1 (after a change of 
notation and correcting typos) reads [16, Eq. (34)] 

N 1 

Z„(l,{mi,m2})« ^ (2Nm 1 m 2 ) 2t+v 

-> I v (^Nm 1 m 2 ) as N -> 00 . (11) 

The final expression Io(AN mi 7712) can be obtained from 
(8) as well, thus verifying our saddle-point result. 

We now comment on the case of odd Nf. With jl = 1, 
assuming m = ml for simplicity, and allowing for v 7^ 
again, we find for the Pfaffians in the integrand of (6) 



Pf 



JY 



Lt _ 



m' 



Pf 



N+v 



m 



m 

L 



(12) 



oc det Ar / 2 (m* 2 l + L^K) ■ det {N+ ^l 2 {m 2 l + LK^) . 

Since K (or L) is an antisymmetric matrix of odd size, at 
least one of its eigenvalues must be zero, and evidently 
the same holds for L^K and LK^ . Thus Z„(l,m) cx 
\m\ 2N m 1 ' , and hence the model with odd Nf seems to 
have an ill-defined chiral limit, with a divergent chiral 
condensate. Further analysis is left to future work. 

While so far we have mainly been concerned with 
the (1 = 1 case, there exists another region of fi which 
is nontrivial, yet analytically tractable: the weak non- 
Hermiticity limit (N > 1 with Nfi 2 fixed) [32], to which 
we now turn for the sake of completeness. First, define 



2) = 



Pt L 



and 971 



m' 
m 



(13) 



In the following discussion Nf is not restricted to be even. 
Since both jl 2 and m are treated as infinitesimal param- 
eters, the saddle point is estimated at ft, — and m = 0, 
i.e., £> = UIU T /4: with U G U(2iV». Taylor expansion 
in jl 2 yields 



Pf 



1 T 



dett/exp [4/t 2 tr(DQ3 T D t$ B) +4tr(D t 2Jl)] 



(14) 



4 



where we used 35^35 = 1/16 and introduced the 2Nf x 

' ^ 1 Isolating the 



2Nf baryon charge matrix 03 _ ( ^ ^ 
U(l) part by U -> Ue l0 / 2N f we find 



d9e lvB / dCf 



SU(2W / ) 



x exp 



rl 



Y/i 2 tr ( UIU T< B 1 \UIU T ) ii % 



2NRetr(e ie/N ?UIU T< m^)} . (15) 



The exponent in this partition function exactly re- 
produces the static part of the chiral Lagrangian [21, 
Eq. (42)] obtained from symmetry principles, with the 
identification \NjjL 2 = V^F 2 ^ 2 and 2Nm = V^F 2 m\ = 
V4\(ipi/j)\M/2Nf, where F is the usual low-energy con- 
stant, m-n- is the pion mass, and in the last equality the 
Gell-Mann-Oakes-Renner relation was used. (A similar 
analysis is given in [20, Eq. (4.15) of 2nd reference] based 
on the one-matrix formulation.) It is intriguing that a 
single chiral RMT, Eq. (4), can describe two extreme 
cases, /i <C Agu( 2 ) and /i s> Agu(2)> that have totally 
distinct patterns of spontaneous symmetry breaking, by 
two different choices of the parameter (p, ~ 0(1/ vN) 
and (1 = 1, respectively) and two different mappings of 
the random-matrix quark masses to the physical quark 
masses (rescaling by |('0V-')l an d A, repectively). 

We expect that this paper will stimulate work in sev- 



eral directions, (i) While we do not doubt that (2) yields 
the correct microscopic spectral correlation functions for 
two-color QCD at high density, for completeness it would 
be useful to prove the equivalence rigorously using the 
partially quenched theory, (ii) We are currently investi- 
gating the calculation of the microscopic spectral correla- 
tion functions using the methods developed in [28]. (iii) 
It would be stimulating to extend our model to two-color 
QCD at intermediate densities, the latter being relevant 
to the phase structure of two-color QCD [20, 33]. The 
issue of continuity between the hadronic phase and the 
BCS superfluid phase is of particular interest, in view of 
the formal similarity of the partition functions at small 
and large /i for degenerate masses [13]. (iv) Since two- 
color QCD with even Nf and pairwise degenerate quark 
masses is free from the sign problem even at /j ^ 0, 
it is in principle possible to test the predictions of our 
model by lattice simulations, (v) Finally, the most im- 
portant extension from a phenomenological viewpoint is 
to the color-superconducting phase of three-color QCD 
[ ]. Work in most of these directions is in progress. 
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